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Abstract 

We study the Schrodinger equation for an off-diagonal tight-binding hamil- 
tonian, as well as the equations of motion for out-of-plane vibrations on the 
separable square Fibonacci quasicrystal. We discuss the nature of the spectra 
and wave functions of the solutions. 

1. Separable Quasicrystal Models 

Much research has been conducted on the properties of excitations in quasicrys- 
tals, ever since the early interest in this question 1 and up to this day. 2 Many exact 
results have been found for 1-dimensional (Id) quasicrystals, yet the properties 
of excitations in 2d and 3d quasicrystals are known to a much lesser degree. Fur- 
thermore, the different models that have been mostly studied — such as Id chains, 
the standard decagonal, octagonal, and icosahedral tilings, or actual structural 
models of real quasicrystals — do not allow one to easily focus on the dimensional 
dependence of the physical properties of quasicrystals, though some interesting 
heuristic arguments have been given by Sire. 3 

One of us has recently studied the geometric properties and calculated the 
diffraction pattern of the square and cubic Fibonacci tilings, suggesting that they 
be used as models for 2d and 3d quasicrystals (with obvious generalization to any 
higher dimension). The advantage of these prototypical models of (i-dimensional 
quasicrystals is that they are separable — certain problems, such as the eigen- 
value problems studied here, can be decomposed into d separate one-dimensional 
problems, yielding straightforward solutions, while allowing directly to focus on 
the effect of dimensionality on the problem being studied. The most apparent 
disadvantage of these models is that they do not occur in "real" quasicrystals, 
yet they should not be dismissed as irrelevant because they can be artificially 
constructed using, for example, conducting nanowires, coupled nanomechanical 
resonators, or photonic quasicrystals. 

The square Fibonacci tiling is constructed by taking two identical grids — each 
consisting of an infinite set of lines whose inter-line spacings follow the well-known 
Fibonacci sequence of short (S) and long (L) distances — and superimposing them 
at a 90° angle, as shown in Fig. 1. This construction can be generalized, of course, 
to any quasiperiodic sequence as well as to higher dimensions. If the original Id 
sequence has inversion symmetry the generated 2d and 3d quasicrystals will have 
square and cubic point group symmetry, respectively. 



Figure 1: A section of the square Fibonacci tiling. 



2. Electrons and (Scalar) Phonons 

Here we take advantage of the separability of the square (and cubic) Fibonacci 
tilings to study two problems: (1) The tight binding hamiltonian with zero on- 
site energy and hopping amplitudes t for vertices connected by long (L) edges, 
1 for vertices connected by short (S) edges, and zero for vertices that arc not 
connected by edges; and (2) The normal modes of out-of-planc vibrations of a 
network of unit masses connected by springs according to the square Fibonacci 
tiling, with spring constants k for long edges, and 1 for short edges. 

The 2-dimensional Schrodinger equation for the connected-neighbor tight- 
binding hamiltonian is given by 

t n+ i*(n + 1, to) + t n ^{n - 1, m) 
+t m+1 $f(n, m + 1) + t m ^(n, m - 1) = EV(n, to), (1) 

where ^(n, m) is the value of a 2d eigenfunction on a vertex labeled by the 
two integers n and to, and E is the corresponding eigenvalue. The hopping 
amplitudes tj are equal to 1 or t according to the Fibonacci tiling as described 
above. The equations of motion for the network of springs is obtained from the 
Schrodinger equation (1) by replacing tj by the spring constant kj, replacing E 
by (fc„+i +k n + k m+ i + k m ) — uj 2 , where u) is the normal frequency, and viewing 
\P(n, to) as the out-of-plane displacement of the (n, to) vertex. The generalization 
of Eq. (1) to three or any higher dimension is obvious. 

With no additional assumptions other than the absence of diagonal hopping 
or diagonal springs this 2-dimensional eigenvalue problem, as well as its higher 
dimensional versions, are all separable. Two-dimensional eigenfunctions can be 
expressed as products 

*ij(n,m) = i)i(n)i)j(m), (2) 

with eigenvalues 

Eij =E t + Ej, or w?j = Luf + lo 2 , (3) 

where Vi( rl ) an d V'i( n ) arc two of the eigenfunctions of the corresponding Id 
eigenvalue equation on the Fibonacci chain, with eigenvalues Ei and Ej or uif 
and u)j, respectively. 
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Separable quasipcriodic hamiltonians have been studied on such models in 
the past, 5 mainly focusing on square or cubic periodic lattices with constant 
hopping amplitudes (or springs) and on-site energies (or masses) that follow the 
same quasiperiodic sequence in all directions. These so-called "diagonal" models 
are separable if one requires the on-site energies (or masses) to be of the form 
V(n, m,...) = V x (n) + V y (m) + . . ., in which case they have solutions as above. 

3. Adding Energies 

When solving the Id tight-binding problem using the standard methods of trans- 
fer matrices and trace maps, 1 one finds for any t that the n th approximant, 
with F n atoms per unit cell, has a spectrum containing F n continuous bands, 
where F n is the n th Fibonacci number. As n increases and the quasicrystal is 
approached, the bands become narrower, and in the limit n — > oo the spectrum 
becomes singular-continuous, containing a Cantor-like set of points whose total 
bandwidth (Lebesgue measure) is zero. Figs. 2(a)-(c) show these bands for the 
first few approximants for three values of the parameter t. One sees that as 
t increases away from the periodic (t = 1) structure the gaps that are formed 
become increasingly wider. 

When adding two such spectra as in Eq. (3) we first note that due to the 
degeneracy of Eij — Eji the maximum number of bands for the n th approximant, 
if there is no overlap at all, is F n (F n + l)/2. The degree of overlap depends on t. 
We find three distinct behaviors which can be described qualitatively as follows: 
(1) For large enough t the gaps in the Id spectrum are sufficiently wide such 
that the 2d spectrum contains well separated bands [Fig. 2(f)]. As n increases 
the total number of bands increases and in the limit n — > oo the 2d spectrum 
becomes singular-continuous with zero total bandwidth. (2) For t close to 1, 
gaps in the Id spectrum are sufficiently small such that the 2d bands greatly 
overlap, forming a single or a few bands even in the n — > oo limit [Fig. 2(d)]. 
The spectrum of the 2d quasicrystal is therefore absolutely continuous. (3) For 
intermediate t a peculiar situation exists (observed also in the diagonal models 5 ) 
where even though the total number of bands increases with n the total integrated 
bandwidth tends to a finite value [Fig. 2(c)]. It is not clear to us at this point 
whether the spectrum of the 2d quasicrystal in this case contains an absolutely 
continuous part or whether it remains singular-continuous as in Id. 

Clearly, the transitions between the different behaviors are pushed to higher 
values of t as the dimension increases and additional Id spectra are added. We 
shall provide elsewhere a detailed analysis of the phase diagram of this system 
as a function of t and dimension. We shall also describe the differences between 
electrons and phonons. 

4. Multiplying Wave Functions 

An intriguing paradox arises in light of the discussion above. Since for any t 
the Id spectrum is singular-continuous, all Id eigenfunctions ipi{ n ) are critical, 
decaying as a power laws from different points no- Any product of two such 
functions, as in Eq. (2), must be critical on the 2d quasicrystal as well, yet we 
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have observed that for t close to 1 the spectrum is continuous, implying that the 
wave functions must be extended. 



(a) Id, t = 1.2 



(b) Id, t = 2.0 



(c) Id, t = 3.0 
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(d) 2d, t = 1.2 
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(f) 2d, t = 3.0 
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Figure 2: Energy spectra of the first few Id and 2d approx- 
imants with t = 1.2, 2.0, and 3.0. In (a)-(c) the number of 
bands for the 3 rd to the 8 th approximants are 3, 5, 8, 13, 21, 
and 34. In (d) there is a single band for all approximants; 
in (e) the number of bands are 3, 3, 9, 13, 19, and 31; and in 
(f) there are 5, 13, 23, 47, 87, and 213 bands. 

We would like to conjecture that the resolution of this paradox stems from the 
high degeneracy of each eigenvalue when bands overlap and merge into one as in 
Fig. 2(d). In general, many different pairs of Id eigenvalues Ei and Ej may add 
up to the same 2d eigenvalue E. Each eigenfunction ^ij(n, m) with eigenvalue 
E is critical, all peaked at different points (no, mo) on the 2d quasicrystal, with 
substantial overlap due to their slow spatial decay. It is plausible that one could 
construct linear combinations of these critical eigenfunctions that are extended 
over the whole infinite quasicrystal. We intend to investigate this conjecture in 



the near future. 

To support our conjecture, we conclude by showing that for any value of t the 
eigenfunctions with energy E = are extended over the 2d quasicrystal. First re- 
call that for any Id eigenfunction ipi{n) with energy Ei, the function (— l) m ipi(m) 
is a Id eigenfunction with energy — Ei. Therefore, ^i{n, m) = (— l) m 7pi(n)'ipi(rn) 
is a 2c? eigenfunction with energy E = 0. Thus, for the n th approximant, the 
energy E — is F„-fold degenerate (a similar situation arises for the labyrinth 
tiling 6 ). Since the Id eigenfunctions form a complete set that spans all func- 
tions on the Fibonacci chain, one can perform a change of basis to an alternative 
complete set <f>j(ri) — J2i c ji^i{ n ) whose members are all extended. The 2d eigen- 
functions 3>j(n, to) = (— l) m <j>j(n)(j>j(m), all with energy E = are extended over 
the whole 2d quasicrystal. 

We should conclude with a cautionary remark that all our conclusions rely 
on the fact that the behavior for finite 2c? approximants survives in the infinite 
limit. This is not obvious until proven rigorously as in Id 
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